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The proximity force approximation (PFA) lias been widely used as a tool to evaluate the Casimir 
force between smooth objects at small distances. In spite of being intuitively easy to grasp, it is 
generally believed to be an uncontrolled approximation. Indeed, its validity has only been tested in 
particular examples, by confronting its predictions with the next to leading order (NTLO) correction 
extracted from numerical or analytical solutions obtained without using the PFA. In this article we 
show that the PFA and its NTLO correction may be derived within a single framework, as the first 
two terms in a derivative expansion. To that effect, we consider the Casimir energy for a vacuum 
scalar field with Dirichlet conditions on a smooth curved surface described by a function ip in front of 
a plane. By regarding the Casimir energy as a functional of tp, we show that the PFA is the leading 
term in a derivative expansion of this functional. We also obtain the general form of corresponding 
NTLO correction, which involves two derivatives of tp. We show, by evaluating this correction term 
for particular geometries, that it properly reproduces the known corrections to PFA obtained from 
exact evaluations of the energy. 



I. INTRODUCTION 

In the last years, there have been important theoretical 
and experimental advances in the analysis of the Casimir 
effect 

Until the recent development of theoretical methods 
that allowed for the exact evaluation of the Casimir en- 
ergy for several geometries, the interaction between dif- 
ferent bodies has been mostly computed using the so 
called proximity force approximation (PFA) [3- This 
approximation, expected to be reliable as long as the in- 
teracting surfaces are smooth, almost parallel, and very 
close, makes use of Casimir's expression for the energy 
per unit area for two parallel plates at a distance a apart. 
For the case of a single masslcss scalar field and Dirichlet 
conditions (the case we deal with in this paper) it is given 
by: 
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The PFA then approximates the interaction between two 
Dirichlet surfaces separated by a gap of spatially varying 
width follows: 



£^PFA = / da Epp{z) , 



(2) 



where S is one of the two surfaces. Quite obviously, this 
formula does not take into account the non-parallelism 
of the surfaces. Moreover, the result may depend on the 
particular surface E chosen to perform the integral. 

As the PFA was believed to be an uncontrolled approx- 
imation, its accuracy has been assessed only in some of 
the particular geometries where it was possible to com- 
pute the Casimir energy numerically or analytically. On 



general grounds, denoting by £ a typical length associ- 
ated to the curvature of one of the surfaces (assumed 
much smaller than the curvature of the second one) and 
by a the minimum distance between surfaces, one expects 
that: 



Ec = EPFA il + j- + 
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where 7 is a constant, whose numerical value fixes the ac- 
curacy of the PFA in each particular geometry (the situa- 
tion could be more complex, since the corrections to PFA 
may contain non-analytic corrections as (^) log(^)). 
One can write similar expressions for geometries that in- 
volve two surfaces of similar curvature. 

In this paper wc explore the following simple idea. 
The Casimir energy can be thought as a functional of 
the shape of the surfaces of the interacting bodies. As 
the PFA should be adequate for almost plane surfaces, 
a derivative expansion |3| of this functional should re- 
produce, to lowest order, the PFA. Moreover, the terms 
involving derivatives of the functions that describe the 
shape of the surfaces should contain the corrections to 
the PFA. We will show that this is indeed the case, and 
that it is possible to find a general formula to compute 
the first corrections to PFA for rather arbitrary surfaces. 

Just to avoid some technical complications, we consider 
a massless scalar field in the presence of a curved surface 
in front of a plane. We will assume that the quantum field 
satisfies Dirichlet boundary conditions on both surfaces. 
Generalizations to other boundary conditions and to the 
electromagnetic field will be analyzed in a forthcoming 
work. 

This paper is organized as follows. In Section |lll we 
describe the model and derive a formal expression for the 
Casimir energy in the geometry described above. Then, 
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in Section lllll wc perform a derivative expansion in the 
expression for the Casimir energy to obtain the main re- 
sult of the paper: a general formula for the interaction 
energy between an arbitrary curved surface and a plane, 
containing up to two derivatives of the function ip that 
describes the curved surface. The leading term of the 
expansion corresponds to the PFA, while the term with 
derivatives is the first non trivial correction. 

In Section IIVI we present some examples: a sphere, a 
cylinder, or a corrugated surface in front of a plane. We 
show, by comparing with existing analytical results, that 
the derivative expansion of the Casimir energy describes 
correctly both the PFA and its first correction for all of 
these geometries. We also compute the derivative ex- 
pansion of the Casimir energy for geometries involving 
parabolic mirrors. Section IVl contains the conclusions of 
our work. 



II. FORMAL EXPRESSION FOR THE VACUUM 
ENERGY 

We shall consider a model consisting of a massless real 
scalar field in 3 -I- 1 dimensions, coupled to two mirrors 
which impose Dirichlet boundary conditions. In our Eu- 
clidean conventions, we use a;o, xi, a;2, ^3 to denote the 
spacetime coordinates, being the imaginary time. 

The mirrors occupy two surfaces, denoted by L and 
R, defined by the equations xa = and 2:3 — 'ip{xi,X2), 
respectively. 

Following the functional approach to the Casimir ef- 
fect, we introduce Z, which may be interpreted as the 
zero temperature limit of a partition function, for the 
scalar field in the presence of the two mirrors. It may be 
written as follows: 



(4) 



where S is the free real scalar field Euclidean action 
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(5) 



while 5l {5r) imposes Dirichlet boundary conditions on 
the L [R) surface. 

Exponentiating the two delta functions by introduc- 
ing two auxiliary fields, Al and Xr, we obtain for Z an 
equivalent expression: 



-S{ip;XL,>.R) 



(6) 



with 



S{^-Xl,Xr)^So{v) (7) 
- i I d*xip{x) [Xl{xi\)S{x3) + XbXxi\)S{x3 - i/'(x||))] 



where we have introduced the notations X|| = (xq, xi, X2) 
and X|| = (xi, X2). 



Integrating out ip, we see that Zq, corresponding to 
the field ip in the absence of boundary conditions factors 
out, while the rest becomes an integral over the auxiliary 
fields: 

2?AL2?Afle"5 /''^^ii E<,,fi-*'-(2;ii)Taf(A3(y||)^ 

(8) 

where a, /? = L, i? and we have introduced the objects: 

TLL(x||,y||) = (x||,0|(-92)-i|y||,o) (9) 

TiH(x||,y||) = (x||,0|(-92)-l|y||,^(y||)) (10) 

T«i(x||,y||) = (x||,V(x||)|(-92)-i|y||,0) (11) 

T^flX^lhyil) = (x||,V(x||)|(-92)-i|y||,V(y||)) (12) 

where we use a "bra-ket" notation to denote matrix ele- 
ments of operators, and is the four-dimensional Lapla- 
cian. Thus, for example. 
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The vacuum energy of the system, -Evac, subtracting the 
zero-point energy of the free field (contained in Zq), is: 

E.,c = lim {^) = T^^TrlogT, (14) 

where T is the extent of the time dimension (or l3~^, in 
the thermal partition function setting), F = — log and 
the trace is meant to act on both discrete and continuous 
indices. 

Note that E'vac still contains 'self-energy' contributions, 
due to the vacuum distortion produced by each mirror, 
even when the other is infinitely far apart. This piece 
(irrelevant to the force between mirrors) shall be sub- 
tracted, in order to obtain a finite Casimir energy, in the 
calculations below. 



III. DERIVATIVE EXPANSION 

We present here a derivation of the first two terms in a 
derivative expansion of the Casimir energy for the system 
defined in the previous section. 

To that end, and for calculational purposes, it is con- 
venient to consider first a simplified situation: we split ip 
into two components. 



(15) 



where a (assumed to be greater than zero), is the spatial 
average of "0, and therefore a constant, whereas rj contains 
the varying piece of ip. The simplified case amounts to 
expanding up to the second order in 77. Since the deriva- 
tives of ip equal the derivatives of 77, to find the terms 
with up to two derivatives of V', it is suflacient to expand 
F up to the second order in 77, keeping up to the second 
order term in an expansion in derivatives: 

F(a,7?) = FW(a) + ^''^\a,T]) + T'^^\a,r]) + ... (16) 
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where the index denotes the order in derivatives. Each 
term will be a certain coefficient times the spatial integral 
over X|| of a local term, depending on a and derivatives 
of rj. 

So far this is a perturbative expansion in rj and its 
derivatives. However, to the same order in derivatives, it 
is quite straightforward to include the corrections which 
arc of the same order in derivatives but of arbitrary order 
in 7]. Indeed, to do this, in the terms obtained in (|16l) . 
one just has to replace a by -0 and also 77 by ijj, before 
performing the spatial integrals. This procedure accounts 
for all the terms of higher order in 77, and the same order 
in derivatives, that contribute to the respective order the 
derivative expansion. Formally, this procedure may be 
represented as follows: 



(17) 



for each term in (|16p . 

Let us calculate the different terms in the derivative 
expansion for F, following this procedure. 

Expanding first the matrix T in powers of 77 



we obtain F = F^") + F^i) + F^^) + 
1. 



(18) 



where 



F(0) = -TrlogT(°) 
F(i) = iTrlog[(T("))-iT« 
r(2) = iTrlog[(TW)-iT(2) 



-^Trlog ^(T(0))-1t(i)(TW)-1t(i) 



(19) 



where, in F^'^ we need to keep up to I derivatives of 77. 

The zeroth-order term is thus simply obtained by re- 
placing first ?/; by a constant, a, and then subtracting the 
contribution corresponding to a — >■ cxj, to get rid of the 
divergent self-energies. This yields, 

F(0)(a) = iTrlog[l-(Tl"2)-iTfi(0-^T^2] (20) 

where the Tjf^ elements are identical to the ones on would 
have for the two fiat parallel mirrors at a distance a apart. 
As mentioned above, we have then to replace a by tp at 
the end of the calculation. After evaluating the trace, we 
obtain: 



T 



log[l 
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(21) 



We then replace a 
Casimir energy. 



(27r)3 

-0 to extract the zeroth order 



d xii 



■ log[l 
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rf^xii 
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(22) 



which equals the PFA approximation to the vacuum en- 
ergy. 

The first order term in the derivative expansion F van- 
ishes identically, while for the second order one we have 
two contributions: 



where, 



p(2) ^ p(2,l) _^ p(2,2) 



r(2a) ^ ^Trlog[(TW)-iT(2) 



and 



r(2,2) = -^Trlog [(T(0))-1t(i)(TW)-1t(i) 



(23) 



(24) 



(25) 



where we have to keep up to two derivatives of 77. 

The form of those terms can be obtained in a quite 
straightforward fashion; indeed, we first note that, in 
Fourier space, and before expanding to second order in 
momentum (derivatives), they have the structure: 



r(2j) 



T 
2" 



<Pk 
(2^)2 



/(2^^)(k)|77(k)|^ 



(26) 



(j — 1, 2), where k = (fci, ^2), f] is the Fourier transform 
of 77, and the /'•^■■'^ kernels are the fco (i.e., static) 
limits of: 



/(^■i)(fc) 

/(2.2)(fc) 



(fp \p\\p + k\ 
(27r)3 1 - e-2|p+fe|a 
(Pp |p||p4- A:|e-2|P+'=l"(l-he-2|pl'i) 

(27r)3 (1 - e-2|p|a)(l _ g-2b+fe|a-) ' 



Besides, we need to subtract an a-independent self-energy 
contribution, obtained by taking a — > 00 in the expres- 
sions above. Putting together the two terms above, and 
subtracting the a — > 00 limit, the total contribution to 
F'^2) adopts the form: 



F 



(2) _ 



T 



(fk 



/(2)(k)|77(k)|^ 



with: 



^ - 'J (27r)3(l 



(2^)2- 

d?p \p\\p + k\ 



g-2|p|a)('g2|p+fc|a _ 



(27) 



(28) 



where we just need to extract its k2 term in a Taylor 
expansion at zero momentum. Namely /'-^'(k) ~ X^^^j 
where 



„ 1 r ^V^^Hfe) - 

^ 2 L 9fc2 ■ 
(Pp 



bl 



lim 



2 r 



\p + k\ 



(27r)3 (1 -e-2|pl'^) k^o dB L(e2|p+''^l" - 1) 
The resulting integral may be exactly calculated, 

^2 



X 



1080 a3 



(29) 
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Thus, 



r(2)(a,,,) 



T 



2 1080 y (27r)2 a3 



2 1080 



(30) 



where, to obtain the second order contribution in deriva- 



tives to the vacuum energy, we need to replace a 
r] ijj, and cancel the T factor, obtaining: 



£;(2) 



r(2)(vj) 



1 



-II 
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T 2 1080 

where the index a runs from 1 to 2. 

Putting together the terms up to second order, 
expression for the energy becomes: 

i?DE = Mac 



^(2) 

vac 
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d^x. 



1 



1 



(31) 



the 



(32) 



This is the main result of this paper. The first term is the 
PFA for the Casimir energy. The second term contains 
the first non-trivial correction to PFA for an arbitrary 
surface. We could have guessed the form of both terms in 
the final formula by using dimensional and symmetry ar- 
guments. The global factor could also be determined by 
considering the particular case of parallel plates. There- 
fore, the calculation presented above, besides confirming 
the general arguments, provides the relative weight be- 
tween both terms, which turns out to be 2/3, regardless 
of the form of the surface. 



IV. EXAMPLES 

We provide here some applications of the general for- 
mula for the Casimir interaction energy. 

A. A corrugated surface in front of a plane 

Let us first consider a corrugated surface in front of a 
plane. For simplicity we assume sinusoidal corrugations 
in the direction of xi 



f 2-KXi 



(33) 



where a is the mean distance to the flat surface, e is the 
amplitude, and A the wavelength of the corrugation. We 
assume a square plane of side L, which is much larger 
than any other length in the problem. 

The derivative expansion for the Casimir energy is 
given by 



E 



DE 




(34) 



In this case, the derivative expansion is an expansion in 
powers of a/A and e/A, i.e. A is largest relevant distance 
in the problem. In order to compare with previous results 
in the literature we will further assume that e <^ a. 
In this limit we obtain 



2t2 r 



E 



DE 
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(35) 



This expression coincides with the small a/ A expansion 
of the result obtained in Ref.Q- Indeed, in that work the 
interaction energy was written as 



E^ 



1440a3 



(36) 



where GtmIx) can be written in terms of Polylogarithm 
functions [5|- One can readily compute the small argu- 
ment expansion of Gtm to obtain 



9 4 9 

TT^ n^x 



480 1080 



(37) 



After inserting this expansion into Eq. p6p, the result 
coincides with the derivative expansion Eq. psp . 



B. A sphere in front of a plane 

We now consider a sphere of radius i? at a distance a 
from a plane. The evaluation of the Casimir energy in 
the electromagnetic case for this configuration has been 
performed in Refs. @, [3|: while the evaluation for scalar 
fields has been previously reported in Ref. See also 
[l3l for asymptotic expansions in the scalar and elec- 
tromagnetic cases near the proximity limit. 

For this geometry, we expect the derivative expansion 
to be adequate in the limit a <^ R. It is worth noting 
that the surface of the sphere cannot be described by a 
single valued function xs = 'il'{xi,X2)- Note that even if 
we consider an hemisphere, the derivatives of will be 
divergent on the equator. For these reasons, the deriva- 
tive expansion will not converge. In spite of this, we will 
see that it still gives quantitative adequate results even 
beyond the lowest order approximation. 

In order to avoid these problems, we will consider only 
the region of the sphere which is closer to the plane. This 
is the usual approach when computing the Casimir en- 
ergy using the PFA. The final result will not depend on 
the part of the sphere considered. Denoting by {p, if) 
the polar coordinates in the {xi,X2) plane the function 
Tp reads 



V'(p) = a + i? 1 - 




(38) 



This function describes an hemisphere when < p < 
R. As mentioned above, the derivative expansion will be 
well defined if we restrict the integrations to the region 
< P < Pm < R- 
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Inserting this expression for tp into the derivative ex- 
pansion for the Casimir energy, one can perform exphc- 
itly the integrations and obtain an analytic expression 
£'De(pm, a, i?). We do not present this rather long ex- 
pression here, but only the leading terms in an expansion 
in powers of a/R, which is given by 

^3 



^(2) 
vac 



R 
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a 

'-R 



and therefore 



R 
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1 a 
SR 



(39) 
(40) 

(41) 



It noteworthy that, up to this order, the result does not 
depend on pu- Moreover, the result is in agreement with 
the asymptotic expansion obtained from the exact for- 
mula for this configuration Q, and with the former nu- 
merical evaluation in [Tl| . 

It is interesting to remark that i^vac includes part of 
the next to leading order corrections. It is correct to keep 
the second term in Eq. (j39p only when the contribution 

(2) 

coming from -Evac is also taken into account. 



C. A cylinder in front of a plane 

Let us now consider a cylinder of radius R and length 
L ^ i? at a distance a from a plane. The Casimir energy 
for this configuration was first evaluated in the PFA in 
Ref ■ [l^ . The exact result was first derived in Ref . [l3j . 
The caveats mentioned in the above subsection also apply 
for this geometry. We will consider the function tjj given 
by 



R 1 



1 



(42) 



with — a^M < xi < xm < in order to cover the part of 
the cylinder which is closer to the plane. The calculation 
is similar to the previous case, and the final result is 



Ebe — — 



I920V2 a^/^ 
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(43) 



Once more, up to this order, the result does not depend 
on xm- Moreover, it is in agreement with the asymp- 
totic expansion obtained from the exact formula for the 
cylinder-plane geometry and numerical findings [l3 - [l7j . 



D. A parabolic cylinder in front of a plane 

We compute here the Casimir interaction energy be- 
tween a parabolic cylinder of length L in front of a plane. 
The surface is defined by the function 



^/'(xi 



■^1 
2J? 



(44) 



with — Xm < Xi < Xm < -R- Once more, we only con- 
sider the portion of the curved surface which is closer to 
the plane (note that the functions defining the cylinder 
Eq. ([^^ and the parabolic cylinder Eq. coincide up to 
first order in xi/ R). The integrations needed to compute 
the derivative expansion of the Casimir energy are very 
simple. Expanding the result in powers oi a/ R we obtain 



1920v^ 



1 



4 a 
9fl 



(45) 



The final answer is independent of xm and the leading 
order coincides with that of the cylinder in front of a 
plane. 



E. A paraboloid in front of a plane 



by 



As a final example we consider a paraboloid, defined 



ibip) = a H 

2R 



with < p < pM < R, in front of a plane. 

The approximation for the vacuum energy reads 



R 



DE 
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(46) 



(47) 



As in all the previous examples, the result does not de- 
pend on the region of integration defined by pu- More- 
over, the leading order is equal to that of the sphere in 
front of a plane, as expected from the fact that the func- 
tions describing both surfaces Eqs. ([55]) and (HH) coincide 
in the region closer to the plane. 



V. CONCLUSIONS 

We have shown that the PFA can be thought of as akin 
to a derivative expansion of the Casimir energy with re- 
spect to the shape of the surfaces. Our main result, given 
in Eq. shows that the lowest order (the "effective 

potential") reproduces the PFA. Moreover, when the first 
non trivial correction containing two derivatives of ip is 
also included, the general formula gives the NTLO cor- 
rection to PFA for a general surface. 

Several remarks are in order: to begin with, at least for 
the surfaces considered in this paper, the PFA becomes 
a well defined and controlled approximation scheme: the 
leading corrections are small when \daip\ 1 or, in other 
words, when the curved surface is almost parallel to the 
plane. Higher order corrections will be negligible when, 
in addition to this condition, the scale of variation of the 
shape of the surface is much larger than the local distance 
between surfaces. It is also clear that the corrections to 
PFA only contains local information about the geometry 
of the surface, and docs not include correlations between 
different points of the surface. 
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Although we apphed our general result to the case of 
a cylinder and a sphere in front of a plane, these ge- 
ometries present additional complications, because they 
cannot be described by a single function -0. Moreover, 
the derivatives of ip diverge when the surface becomes 
perpendicular to the plane, and therefore it is clear that 
the derivative expansion will not converge. In spite of 
this, it is remarkable that Eq. ([5^ describes the inter- 
action energy for these configurations including the first 
non trivial correction to PFA. Strictly speaking, for these 
geometries we are computing the interaction energy be- 
tween a plane and a large curved surface which, in the 
region closest to the plane, has a cylindrical or spherical 
shape. 

We expect the main idea presented in this paper to be 
gencralizable in several directions, as for instance for a 
scalar field satisfying Neumann or Robin boundary con- 
ditions, and also to the electromagnetic field satisfying 
perfect conductor boundary conditions on the surfaces. 



In all these cases, we expect the derivative expansion to 
be of the form 

= " llio / '^'''11 ^ + , (48) 

where the constants /?; will depend on the kind of fields 
and boundary conditions considered. 

Other interesting generalizations would be to consider 
two curved surfaces, and the case of imperfect boundary 
conditions. Moreover, as the applications of the PFA 
are not restricted to the Casimir energy, the derivative 
expansion could also be useful to compute gravitational 
, electrostatic [l^ or even nuclear forces [1^ . 
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